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ABSTRACT 

Using the Green Bank Telescope (GBT), we have obtained accurate measurements of the ^"^N/^^N 
isotopic ratio in ammonia in two nearby cold, dense molecular clouds, Barnard 1 and NGC 1333. The 
i4]\j^i5]\j ratio in Barnard 1, 334 ± 50 (3(t), is particularly well constrained and falls in between the 
local interstellar medium/proto-solar value of '^450 and the terrestrial atmospheric value of 272. The 
NGC 1333 measurement is consistent with the Barnard 1 result, but has a larger uncertainty. We do 
not see evidence for the very high ^'"'N enhancements seen in cometary CN. Sensitive observations of 
a larger, carefully selected sample of prestellar cores with varying temperatures and gas densities can 
significantly improve our understanding of the nitrogen fractionation in the local interstellar medium 
and its relation to the isotopic ratios measured in various solar system reservoirs. 
Subject headings: astrochemistry — ISM: abundances — ISM: individual (Barnard 1, NGC 1333) — 
ISM: molecules — molecular processes — radio lines: ISM 



L INTRODUCTION 

Molecular isotopic ratios are an important tool in the 
study of the origin of interstellar and solar system materi- 
als. After hydrogen, nitrogen displays the largest isotopic 
variations in the solar system, typically explained by mix- 
ing of various proto-solar or pre-solar reservoirs. Earth, 
Mars interior, Venus and most primitive meteorites have 
nitrogen isotopic ratios within 5% of the terrestrial atmo- 
spheric value, "N/i^N = 272 (see lMartv et al. II2OO90 . 
However, the proto-solar nebula was poorer in ^^N, with 
^^N/^^N ~ 450, as evidenced by infrared and in-situ 
measurem ents i n the Jupiter atmosphere ( 530liyn, ISO, 
Fouch et et al. I 120001: 4 35 ± 57, Galileo. lOwen et al. I 
200ir4 48 ± 62. Cassini. lAbbas et al. I |20M; 450 ±106^ 
Cassini. lFouchet et al. \\2004) and recent solar wind mea- 
surements (442 ± 131, 2cr, Genesis, iMartv et al~ll2009f ). 
The proto-solar ^*N/^^N ratio is in agreement with 
the local interstellar medium (ISM) value (450 ± 22, 
IWilson fc Rood II1994I: or 414 ± 32 at the birth place of 
the Sun. lWielen fc WilsoTTIfTool . 

Very low nitrogen isotopic ratios, ^"'N/^^N = 148 ± 
6, have been mea sured in C N in a large sample of 
comets (j^'Ianfroid et al. I l2009f ) . Similar isotopic ratios 
have been recently derived by iBockelee-Morvan et al. I 
(|2008[) in comet 17P/Holmes in CN and HCN, the pre- 
sumed parent of the CN radical in cometary atmospheres 
("N/i^N = 165 ± 40 and 139 ± 26, respectively). Al- 
though higher ratios in HCN were initially reported in 
comet Hale-Bopp, Bockelee-Morvan's re-analysis of the 
archival data suggests that the -'^"'N/^'^N ratio in HCN in 
this comet is also consistent with the CN value, given 
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the measurement uncertainties. The exact chemical net- 
works explaining the HCN isotopic anomaly in comets 
have yet to be proposed. However, one clear conclu- 
sion is that HCN has never isotopically equilibrated with 
the nebular N2 at the later phases of the evolu tion of 
the solar system (|Bockelee-Morvan et al. II2008D . Low 
nitrogen isotopic ratios are also found in some mete- 
orites, interplanetary dust particles (IDPs) and samples 
of material returned from comet 81P/Wild 2 by the Star- 
dust spacecraft; am ine groups in pa.rticular often show 
low ^ " ^N/^^N ratios fEngel fc Macko Ill997t iKeller et al. 1 
[200I iFloss et al. 1 12006: McKeeganet al. il2006D . 

In the ISM, nitrogen bearing molecules are important 
tracers of cold, dense g they do not freeze-out 

onto grain mantles at densities below a few xlO^ cm~^, 
contrary to carbon and oxygen species, such as CO 
or CS. In cold, dense, CO depleted ISM regions, high 
^^N enhancements have been suggested to be present 
in gas-phase molecules, in pa rticular in ammonia (see 
iRodgers fc Charnlev I l2008aU bl. and references therein). 
These isotopic enhancements could be preserved as spa- 
tial heterogeneity in ammonia ice mantles, as monolayers 
deposited at different times would have different isotopic 
compositions and may in turn explain the ^''N enhance- 
ments seen in primitive solar system refractory organics, 
if they are synthesized from ammonia. 

Measurements of nitrogen isotopic ratios in ammonia 
in cold molecular clouds are of great interest, as they pro- 
vide quantitative observational constraints for the chem- 
ical fractionation models. However the data a vailable are 
scarce. In a recent study, iGerin et al. I (|2009[ ) detected a 
doubly-fractionated ammonia isotopologue, ^^NH2D, in 
several ISM regions. They derive ^^N/^^N isotopic ratios 
~350-850, similar to the proto-solar value of 450. How- 
ever, the ^^NH2D lines are weak and the corresponding 
measurement uncertainties are large. To better constrain 
the nitrogen isotopic ratios in the solar neighborhood, we 
have observed the -'^"'NHa and -'^^NHs inversion lines using 
the NRAO 100 m Robert C. Byrd Green Bank Telescope. 
We present here results for two nearby clouds, Barnard 1 
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Fig. 1. — Spectra of the ammonia inversion lines in Barnard 1- 
b (aj2000 = 03''33'"20f8, Sj2ooo = +31°07'34"). Gray lines 
show hypcrfine structure fits — a single velocity component for the 
^"'NHa (1,1), (2,2), and l^NHa (1,1) transitions and two compo- 
nents for the ^''NHa (3,3) transition, as described in the text. The 
^^NHa (2,2) spectrum, shown with a Gaussian line fit, is formally a 
3.5(7 upper limits. The bottom-right panels shows the LTE rotation 
diagram, as described in the text. LHS stands for the left-hand-side 
in eq. (2) of Blake et al. (1987). 

and NGC 1333, and discuss prospects for future studies 
of nitrogen isotopic ratios using the GBT. 

2. OBSERVATIONS 

Observations of ammonia inversion lines presented here 
were carried out in 2009 February using the GBT. Con- 
ditions were calm under a lightly overcast sky resulting 
in a stable system temperature of '^48-50 K. The K- 
band receiver feed/amplifier set covering 22-26.5 GHz 
was employed. The GBT spectrometer was configured 
in its eight 50 MHz intermediate frequency (IF) band- 
width, nine-level mode, which allows simultaneous cover- 
age of four 50 MHz frequency bands in two polarizations, 
through the use of offset oscillators in the IF. This mode 
gives 3.051 kHz channel separation. The IF No. was 
centered between the (1,1) and (2,2) lines of ^''NHa, the 
(3,3) line was centered in IF No. 1 and the (2,1) line in 
IF No. 3; IF No. 2 was tuned midway between the (1,1) 
and (2,2) lines of ^^NHg. 

The antenna temperature calibration was achieved 
by injecting a calibrated, broadband signal from a 
noise diode. The d ata were placed on the scale 
(jUhch fc Haas I 119761 ) by multiplying by a factor of 
e'^^/rji, where r is the zenith opacity (^-^0.065 at 23.5 GHz 
and 0.075 at 22.5 GHz, from an atmospheric model), A 
is the air mass of the observations, and rji is the rear 
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Fig. 2. — Spectra of the ammonia inversion lines in NGC 1333 

(aj2000 = 03h29-"ll?6, 5j2ooo = +31°13'26"). The i^nHs (2,1) 
spectrum is formally a 3.7cr upper limit. The bottom-right panel 
shows the LTE rotation diagram, as described in the text. Lines 
as in Fig. 1. 

spillover efficiency (^0.99 for the unblocked GBT aper- 
ture). The source 0336+3218 was used for pointing and 
absolute calibration at the beginning of the observation 
set. The absolute calibration accuracy is conservatively 
estimated at 20%. Relative accuracy among lines in 
spectral bandpasses measured at the same time is much 
better, typically dominated by the random noise of the 
spectra. The conventional beam efficiency of the GBT 
at 20 GHz is r]mb ^ 0.79 and the FWHM beam size is 
~33". Data were taken in the frequency switching mode, 
with a frequency displacement of 6 MHz. 

Automatically updated dynamic pointing and focusing 
corrections were employed based on real-time tempera- 
ture measurements and a thermal model of the GBT; 
zero points were adjusted typically every 2 hours or less, 
using 0336-1-3218 as a calibrator. The two polarization 
outputs from the spectrometer were averaged in the fi- 
nal data reduction process to improve the signal-to-noise 
ratio. The total integration times for Barnard 1 and 
NGC 1333 are 174 and 48 min, respectively. 

3. RESULTS 

Figure [T] shows spectra of the ^"^NHs and ^^NHa in- 
version lines toward Barnard 1. The satellite hyper- 
fine lines arc detected with high signal-to-noise ratios 
for the (1,1) and (2,2) transitions of ^^NHg. Resuhs of 
hypcrfine structure (HFS) fits, obtained using the IRAM 
CLASS software package, are shown in Tablc[TJ The (1,1) 
and (2,2) ^"'NHg lines have the same line width of 0.80 
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TABLE 1 

HFS Fit Results and Ammonia Column Densities. 



Transition 




Vlsr 


AV 


T 


N 1 




(K km s-i) 


(km s-i) 


(km s-i) 




(cm^^) 


Barnard 1-b 


"NHa (1,1) 


37.7 ±0.006 


6.7 


0.80 


4.0 


2.8 X IQiS 


"NHs (2,2) 


1.90 ±0.005 


6.5 


0.80 


0.3 


2.8 X IQis 


"NH3 (3,3) 


0.076 ± 0.029 


6.5 


O.SO*" 




8.0 X 10^^ 




0.26 ±0.11 


7.1 


3.5 






"NHs (2,1) 


0.012 ± 0.004 










15NH3 (1,1)'= 


0.106 ±0.004 


6.6 


0.63 




8.4 X 10^2 


i^NHs (2,2) 


0.013 ±0.004 










NGC 1333 


"NH3 (1,1) 


16.4 ±0.015 


7.5 


1.2 


1.2 


8.9 X 10" 


"NHs (2,2) 


2.22 ±0.07 


7.3 


1.3 


1.1 


9.9 X 10" 


"NH3 (3,3) 


0.37 ±0.06 


6.8 


1.3'' 




2.6 X 10^5 




0.61 ±0.04 


4.7 


9.5 






"NH3 (2,1) 


0.023 ±0.007 










i^NHs (1,1)^= 


0.045 ±0.007 


7.2 


0.85 




2.6 X 10^2 


i^NHg (2,2) 


0.012 ±0.008 











Notes: Upper limit for the ^^NHs (2,1) and ^^NHa (2,2) transitions included for completeness. Molecular column densities are computed 
assuming LTE with temperatures of 10.9 and 14.7 K for Barnard 1 and NGC 1333, respectively, as described in the text. "Uncertainties are 
from the HFS fit for the ^^NHs (l,l)-(3,3) lines and computed from the rms in the spectra for the remaining weak lines and upper limits. 
''Fixed parameter. "^The hyperfine structure of ^^NHa has been recently revised bv lBethlem et al. I I I200SI ) who give the lev el energies with 
unprecedented accuracy. We have calculated the corresponding transition strengths using the theory of IKukolichl ||T967| ') and assummg 
half-integer values of the quantum numbers of the levels involved. 



kms~^, while the (3,3) spectrum is much broader, indi- 
cating the presence of a high-velocity component origi- 
nating in warmer gas. Consequently, we fitted the (3,3) 
line with two velocity components, fixing the line width 
of the first component at the value obtained for the (1,1) 
and (2,2) lines. All parameters for the second component 
were allowed to vary in the fit. The ^^NHs (1,1) line is 
slightly narrower than the ^^NHa (1,1) and (2,2) lines, as 
expected for an optically thin line. Figure [2] shows the 
NGC 1333 spectra. The lines are weaker and broader 
than in Barnard 1. The (3,3) transition shows strong 
blueshifted emission from the molecular outflow driven 
by IRAS 4A. The line width derived from the HFS flt to 
the ^^NHa data is significantly smaller than the values 
obtained for ^^NHs. However, the signal to noise ratio 
in the ^^NHa spectrum is relatively low, formally 6.2(t. 

We use the rotation diag ram technique (see, e.g., 
eq. 2 of IBlake et al. I Il987f ) to derive the excitation 
temperatures and molecular column densities. We use 
the spectroscopic constants and partition functions for 
^''NHa an d ^^NHs from^ the JPL Mol ecular Spectroscopy 
database ( http: / /spec.jpl.nasa.gov/[ ) . The resulting ro- 
tation diagrams are shown in l^ igureslTland O bottom- 
right panels. The least-squares fit to the line intensities 
of the (l,l)-(3,3) data (narrow component) in Barnard 1 
gives a rotation temperature of 12.4 K in (dotted line 
in Fig. [1]). However, since the (3,3) spectrum is con- 
taminated by the high- velocity emission, only the (1,1) 
transition is detected for ^^NHa, and given our interest in 
the cold gas component on the line of sight, a more con- 
sistent approach is to use only the (1,1) and (2,2) ^"'NHa 
data for the temperature determination. This approach 
also mitigates p ossible effects of the no n-LTE ortho-to- 
para ratio (e.g.. lUmemoto et al. Ill999t) and results in a 
slightly lower rotation temperature of 10.9 ± 0.2 K (as- 



suming 2.7% uncertainties for the integrated line inten- 
sities, as discussed in § |4l dashed line in Fig. [1] bottom- 
right panel). This value is subsequently used to derive 
the ^^NHa and ^^NHa column densities from the (1,1) 
integrated line intensities (Table [1]). Given the low rota- 
tion temperature, we include the correction for the 2.7 K 
cosmic background radiation in the column density com- 
putations. This correction has no effect on the tempera- 
ture determination, since all the lines observed have very 
close rest frequencies. The NH3 temperature and column 
density in Barnard 1 we derive here are i n good agree- 
ment with the results of iBachiller et al. I (|1990[ ) based 
on observations with the Effelsberg telescope (12 K and 
2.5 X IQis cm-2). 

For NGC 1333, we derive a rotational temperature of 
17.6 K using the (l,l)-(3,3) lines and 14.7 ±0.5 K using 
only the (1,1) and (2,2) transitions (dotted and dashed 
lines, respectively in Fig. [21 bottom-right panel). We use 
the latter value for the column density determination. 

4. DISCUSSION 

The uncertainty of the derived column densities is dif- 
ficult to estimate quantitatively. As discussed in § [51 all 
lines are observed simultaneously and the relative cali- 
bration is largely determined by the signal-to-noise ratio 
in the spectra — 24 and 6.2ct for ^^NHa (1,1) in Barnard 1 
and NGC 1333, respectively. The "NHg/^^NHs column 
density ratio, when computed from the observed inten- 
sities of the (1,1) transitions of the two isotopologues, is 
insensitive to the exact value of the rotation tempera- 
ture used in the calculations, as long as the temperature 
is the same for the two isotopologues. The remaining 
source of uncertainty is the determination of the ^"^NHa 
optical depth. The formal uncertainties for the opacity 
corrected line intensities given by the HFS fits are negligi- 
ble (see Table [J). In order to independently estimate the 
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uncertainty of the ^^NH3(1,1) integrated line intensities 
obtained from the HFS fit, we divide the full Barnard 1 
data set into three subsets, 58 min integration time each, 
and fit each subset independently. From the three mea- 
surements we derive a 2.7% uncertainty for the mean 
which we use as an estimate of the uncertainty of our fi- 
nal "'^^NH3(1,1) measurement. Adding in quadrature the 
corresponding statistical uncertainty for ^^NHg of 4.2%, 
we derive a Scr uncertainty for the ^"'NHaZ-'^^NHa ratio of 
15.0% in Barnard 1. Assuming a 5.2% uncertainty for the 
^"'NHa (1,1) measurement in the shorter NGC 1333 in- 
tegration (the Barnard 1 value scaled by the square root 
of the integration time ratio) , and adding in quadrature 
a 16% uncertainty for the ^^NHa measurement, we de- 
rive a 50% uncertainty (3cr) of the i^^NHa/i^NHg ratio 
in this source. Our final estimates of the -'^"'NHs/^'^NHa 
abundance ratio are thus 334 ±50 and 344 ±173 (3cr) for 
Barnard 1 and NGC 1333, respectively. 

The ^^N/^^N fractionation ratios in Barnard 1 and 
NGC 133 3 derived from deu terated ammonia measure- 
ments bv lGerin et aFI (12001) are 470t}^j] and 360t^|[!, 
respectively. The NH3 and NH2D results are thus consis- 
tent, given the NH2D measurement uncertainties. Owing 
to the relatively short integration and low SNR in the 
^^NHs spectrum, our result for NGC 1333 is not very 
constraining. However, our Barnard 1 data are of excel- 
lent quality. The derived ^*N/^^N ratio falls in-between 
the local ISM/proto-solar value (~450, as discussed in 
§ 1) and the terrestrial atmospheric ratio of 272. We do 
not see evidence for the very high ^^N enhancements, as 
suggested by the cometary CN measurements (^-^150). 

A comparison of our fractionation measurements 
with the nitrogen "super fractionation" models of 
iRodgers fc Charnlev I (l2008bl) is of in terest. In their Fig- 
CharnlevT (l2008bO present results of 
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time-dependent models of the evolution of the ^^N/^'^N 
enhancement ratios in key species, including ammonia. 
For a temperature of 10 K, the gas-phase ^'^N fraction- 
ation in ammonia can reach a maximum enhancement 
of a factor ~7-10, much higher than the value we de- 
rive for Barnard 1. However, the enhancement factor 
decreases with time to a factor of a few below the initial 
isotopic ratio assumed in the model. The time at which 
the maximum fractionation is reached depends strongly 
on the initial fraction of nitrogen in molecular form, but 
is typically of order 1 Myr. The enhancements of the 



magnitude that we measure can be present very early 
on, or a few tenths of a Myr after the peak enhancement 
time. 

The IRodgers fc Charnlev I ()2008bl ) model rcsuhs de- 
pend strongly on the gas temperature — at 7 K, sig- 
nificantly higher enhancements are found in the gas- 
phase ammonia fractionation, as compared to 10 K. 
The sources we have studied are s omewhat warmer than 
the mo del clouds considered by IRodgers fc Charnlev I 
(|2008b( ) . It would thus be extremely interesting to extend 
our measurements to colder clouds, such as LDN 1544 
or LDN 134N, with central temperatures of order 7 K. 
Our present observations are exploratory in nature, but 
show that accurate measurements of nitrogen fractiona- 
tion in ammonia in cold ISM sources are clearly feasible 
using the GBT. A focused observing program of a larger, 
carefully selected sample of prestellar cores with varying 
temperatures and gas densities would allow significant 
progress in our understanding of nitrogen fractionation 
in local ISM and its possible implications for the existing 
measurements in various solar system reservoirs. 

The nitrogen fractionation enhancement levels derived 
here arc well reproduced by the gas-phase models of 
Roucff ct al. (see, e.g.. Fig. 2 of Gerin et al. 2009). 
For an H2 density of 10^ cm~'^ and a 10 K temperature, 
the gas-phase models predict ^^NHs/^^NHa = 389 and 
i^NHsD/i^NHaD = 372. The models of Roueff et al. 
are steady-state models that focus on the gas-phase pro- 
cesses and approximate the adsorption of the species on 
grain surfaces by reducing the elemental abundances in 
the gas, as compared to the dynamical models of Rodgers 
fc Charnlcy that explicitely take into account the mantle 
formation and gas-grain interactions. 
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